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We present a theoretical analysis of the electron spin decoherence in single Nitrogen- Vacancy de- 
fects in ultra-pure diamond. The electron spin decoherence is due to the interactions with Carbon-13 
nuclear spins in the diamond lattice. Our approach takes advantage of the low concentration (1.1%) 
of Carbon-13 and their random distribution in the diamond lattice by an algorithmic aggregation 
of spins into small, strongly interacting groups. By making use of this disjoint cluster approach, we 
demonstrate a possibility of non-trival dynamics of the electron spin that can not be described by a 
single time constant. This dependance is caused by a strong coupling between the electron and few 
nuclei and results, in particular, in a substantial echo signal even at microsecond time scales. Our 
results are in good agreement with recent experimental observations. 



I. INTRODUCTION 

Isolated spins in solid-state systems are currently 
being explored as candidates for good quantum bits, 
with applications to quantum computational 2 ^, quan- 
tum communication^ and magnetic sensing 5 -! 6 -! 7 -. The 
Nitrogen-vacancy (NV) center in diamond is one such 
isolated spin system. It can be prepared and detected 
using optical fields, and microwave radiation can be used 
to rotate the spin 8 -! 9 -. Recent experiments have conclu- 
sively demonstrated that in ultra pure-diamond the elec- 
tron spin coherence lifetime is limited by its hypcrfinc 
interactions with the natural 1.1% abundance Carbon-13 
in the diamond crysta l 10 ! 11 . Thus, developing a detailed 
understanding of the decoherence properties of such an 
isolated spin in a dilute spin bath is a challenging prob- 
lem of immediate practical interest. This combined sys- 
tem of electron spin coupled to many nuclear spins has a 
rich and complex dynamics associated with many-body 
effects. 

The decay of electronic spin coherence due to inter- 
actions with surrounding nuclei has been a subject of 
a number of theoretical studie o 12 ' 13 . Various mean-field 
and many-body approaches have been used to address 
this proble m 14 i 15 i 16 i 17 ' 18 i 19 . In this paper, we investi gate 
a variation of the cluster expansion, developed in Ref. Ilq . 
Our approach takes advantage of the natural grouping 
statistics for randomly located, dilute impurities, which 
leads to the formation of small, disjoint clusters of spins 
which interact strongly within themselves and with the 
central spin, but not with other such clusters. This sug- 
gests a natural hierarchy of interaction scales of the sys- 
tem, and allows for a well-defined approximation that can 
be seen as an extension of ideas developed in the study 
of tensor networks 20 . We develop an algorithm for find- 
ing clusters given a set of locations and interactions, and 
find that for dilute systems convergence as a function of 
the cluster size (number of spins in a given cluster) is 
very rapid. We then apply this technique to the partic- 
ular problem of the decay of spin-echo for a single NV 



center, and find good qualitative and quantitative agree- 
ment with experiments. In particular, we demonstrate 
a possibility of non-trival dynamics of the electron spin 
that can not be described by a single time constant. This 
dependance is caused by a strong coupling between the 
electron and few nuclei and results in a substantial spin- 
echo signal even at microseconds time scale. 



II. METHODS 

In this section, we introduce the properties of the elec- 
tron spins in a NV center and the nuclear spins in its 
environment. Then, we address the many body problem 
involved in the evaluation of spin-echo signals. 



A. Spin hamiltonian 

The negatively charged NV center ([N-V] - ) has trig- 
onal C$ v symmetry and 3 A2 ground stated with total 
electronic spin S = l 22 . Spin-spin interaction leads to a 
zero-field splitting, A = 2.87 GHz, between the m s = 
and m s = ±1 manifolds, where the quantization axis is 
along the NV-axis. This spin triplet interacts via hyper- 
fine interaction with a spin bath composed of the adjacent 
Nitrogen-14 and the naturally occuring 1.1% Carbon-13 
which is randomly distributed in the diamond lattice. 

In the presence of an external magnetic field, the dy- 
namics is governed by the following hamiltonian, 

H = A5 2 - 7e B 2 5 z -^ 7w B.g„(|5 2 |)-I„ 

n 

+ Y / S z A n -I n + Y / SA n (\S z \)-I n 

n n 

+ ^I„-C nm (|S 2 |)-I m . (1) 

n>m 

The relatively large zero-field splitting A (first term in 
Eq. iJTJl) does not allow the electron spin to flip and 



thus we can make the so called secular approximation, 
removing all terms which allow direct electronic spin flips. 
Non-secular terms have been included up to second order 
in perturbation theory, leading to the \S Z \ dependence of 
other terms in the Hamiltonian. The second and third 
terms are, respectively, the Zeeman interactions for the 
electron and the nuclei, the fourth term is the hyperfine 
interaction between the electron and each nucleus, the 
fifth term is an effective crystal-field splitting felt by the 
nuclear spins, and the last term is the dipolar interaction 
among nuclei. The specific terms for this hamiltonian are 
discussed in the Appendix. 

For the NV center case, the nuclear g-tensor, g„, can 
be anisotropic and vary dramatically from nucleus to 
nucleus-^. This leads to a non-trivial dynamics between 
the electron and an individual nucleus (electron-nuclear 
dynamics) , and motivates a new approach for the case of 
a dilute bath of spins described below. In addition, the 
interaction between nuclei is enhanced by the presence 
of the electron of the NV center. The resulting effective 
interaction strength can exceed several times the bare 
dipolar interaction between nuclei^. 



B. Disjoint cluster method 

The large zero-field splitting A sets the quantiza- 
tion axis (called NV-axis) and allows us to neglect elec- 
tron spin flips due to interactions with nuclei. There- 
fore, we can reduce the Hilbert space of the system by 
projecting hamiltonian (|A3[) onto each of the electron 
spin states. We can write the projected hamiltonian, 
P ms HP ms (where P ms = \m s ){m s \), as 

H ms = ^^)-I n +^I n -Cl^)-I ro +A|m 5 |- 7e -B z m SI 

n nm 

(2) 

where m s denotes the electron spin state, i~2„ is the effec- 
tive Larmor vector for nucleus n and C nrn is the effective 
coupling between nuclei n and m. In equation @, we in- 
clude the zero-field splitting and the Zeeman interaction. 
These terms provide just static fields whose effect is can- 
celed by spin echo. In this way, we can write the evolu- 
tion of the bath as U„ ls (r) = T {exp (/J" H m . (*') dt') }. 
An exact expression for U ms can be found by ignoring 
the intra-bath interactions (C„ m = 0). However, for an 
interacting bath (with arbitrary C nm ), solving U ma for 
a large number of nuclei N is a formidable task since 
it requires describing dynamics within a 2 N dimensional 
Hilbert space. Therefore, some degree of approximation 
is needed. 

The spin bath considered here is composed of randomly 
distributed spins, and not all pair interactions among 
nuclei are equally important. Specifically, interactions 
decay with a characteristic law 1/R^ nn , where R nm is 
the distance between nuclei n and m. As a result, we 
can break the big problem into smaller ones by group- 
ing those nuclei that strongly interact with each other 
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FIG. 1: Illustration of the method. Spins that strongly inter- 
act can be grouped together and treated as isolated systems. 
Interactions that joint different groups can be incorporated 
as a perturbation. 

into disjoint sets. Our procedure is illustrated in figure 
Q] We denote the k-th group of nuclei as Cg , where the 
subindcx g indicates that each group has no more than g 
nuclei. Interactions inside each group (intra-group inter- 
actions) are expected to be much larger than interactions 
among groups (inter-group interactions). Our approxi- 
mation method will rely upon neglecting the latter. 

Formally, we start by separating intra-group interac- 
tions and inter-group interactions. We define the opera- 
tor Hb = Y2k H(Cg) which contains all electron- nuclear 
interactions (first term in Eq. ((2|)) plus all interactions 
between bath spins within the same group Cg. Simi- 
larly, operator Ha = H(C g )(= H — Hb) contains all in- 
teractions between bath spins in different groups. As 
a first approximation, we can neglect the inter-group in- 
teractions but keep the intra-group interactions. The ap- 
proximation can be understood by means of the Trotter 

24 

expansion 

exp (HaI + Hst) = lim (exp (HAt/n) exp (i^ei/n))™ . 

n — >oo 

Since Hb contains groups of terms that are disconnected 
from each other, [H(Cg), H(Cg )] = and we can write 
the evolution operator as 

u 'V=£z. (^>3n^D) > ( 3 ) 

where U(C g , — ) is the evolution operator due to hamilto- 
nian H{Cg) and so on. 

To the zcroth order we neglect all terms in Ha since 
H(C g ) contains interactions among nuclei that interact 
weakly. Thus, we set U(C g , — ) to the identity and sim- 
plify Eq. © to 

Ug(r)^Y[u(C k g ,r). (4) 

k 

This approximation requires independent calculation of 
propagators for each group g, which corresponds to N/g 



3 



independent calculations of 2 9 x 2 9 matrices, exponen- 
tially less difficult than the original problem of direct cal- 
culation of the 2^ x 2 N dimensional matrix. We remark 
that including the effect of U(C g , — ) in the Trotterization 
can be done by using a tree tensor network ansatz wave 
function^ where the number of complex coefficients to 
describe the wave function is 0(N l ° 9 ^) instead of 2 N . 



III. ELECTRON SPIN-ECHO 

Electron spin-echo removes static magnetic shifts 
caused by a spin bath, allowing to measure the dynami- 
cal changes of the bath. Assuming that an initial state, 
If) — (|0) + |1})/V2) is prepared, the probability of re- 
covering the same state after a time 2t is 



P = Tt(P v Ut(t) P U^(t)), 



(5) 



where P v — \ip)(cp\ is the projector operator to the ini- 
tial state, p = \ip)((p\ (g> p n is the density matrix of the 
total system, p n is the density matrix of the spin bath, 
Ut (t) = U (t) RttU (t) is the total evolution of the sys- 
tem where U is the evolution operator under hamilto- 
nian ([1]) and is a 7r-pulse acting on the subspace 
m s = {0, 1} of the electron spin manifold. Probability 
([5]) can also be written as p = (1 + S (t))/2, where 



S (r) = Tr (p n ul (r) U\ (r) U (r) U, (r) 



(6) 



is known as the pseudo spin and |5(r)| = is the long- 
time (completely decohered) signal. In the high tem- 
perature limit, the density matrix of the nuclei can be 
approximated by p n ~ 1® N /2 N where N is the num- 
ber of nuclei. The generalization of this relation for 
different sublevels of the triplet state is straighforward, 

S (r) = Tr [p n Ul (r) u\ (r) U a (r) Up (r)) , where a = 1 

and (3 = — 1, for example. In what follows, we analyze 
the effect of an interacting bath on Eqn. ©. 



When the electron spin is in its m s = state and the 
external magnetic field points parallel to the NV-axis, 
the Larmor frequency ft„ is set by the external mag- 
netic field, and the nuclei precess with the same frequency 
The total pseudo spin is 1 at times ft^V — 2m7r 
with m integer. When the electron is in its m s = 1 
state, the Larmor vector fin has a contact and dipolar 
contribution from the hyperfine interaction A„ that may 
point in different directions depending on the position of 
the nucleus. As a consequence, when interactions from 
all nuclei are considered, these electron-nuclear dynamics 
makes the total pseudo-spin relation collapse and revive. 
However, it does not show any decay of the revival peaks. 

We point out that when the transverse (perpendicular 
to the NV-axis) magnetic field is non-zero, nuclei near 
the center experience an enhancement in their g-factors 
leading to a position-dependent Larmor frequency ftl°^ 
(see the Appendix). This will result in an effective decay 
of the signal since the electron state will not be refocused 
at the same time for all nuclei. 



B. Interacting bath: an example 

When the intra-bath interactions are considered, the 
spin-echo signal can show decay in addition to the 
electron-nuclear dynamics. As an illustrative and sim- 
ple example, consider a pair of nuclei with their Larmor 
vectors pointing in the same direction regardless the elec- 
tron spin state (in this case there is no electron-nuclear 
dynamics and the non-interacting pseudo-spin relation 
for two nuclei is S nm = S n S m = 1 (see Eq. 0)). When 
the interaction between nuclei is included, the pseudo- 
spin relation can be worked out exactly, 
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(8) 



A. Non-interacting bath 

To understand the effect of an interacting bath we will 
first analyze the non-interacting case, which displays the 
phenomenon of electron spin-echo envelope modulation 
due to electron spin- nuclear spin entanglement. This is 
completely neglecting interactions among nuclei, C nm — 
0. In this regime, the evolution operator is factored out 
for each nucleus and the pseudo-spin is the product of 
all single pseudo-spin relations. In the high temperature 
limit, p n = 1/2, we obtain the exact expression^ 



where (u^O = (Aft^/2) z + fe) 
ft™ 5 and is the strength of the dipolar interaction 



Aft™ 3 — ft m = — 

'-"'nm — iL n 



C nn\{^n+Im- 



I n. — II 



m+ 



i ) between nuclei n and 



to. The two frequencies involved in ([5]), w nm and w^ m , 
are not necessarily the same for different pairs of nu- 
clei. They depend on the relative position between nu- 
clei and the relative position of each nucleus to the NV 
center. Therefore, when all pair interactions are included 
the pseudo-spin relation decays. In the following section 
we present an approach to incorporates not only this two 
body interaction but also n-body interactions with n < 6. 



s T (t) = n 5 «M=n( 1 - 2 



ft (0) T 

x sin sin 



ftl°) X ft« 



ft ( V 



(7) 



IV. THE DISJOINT CLUSTER APPROACH 

The many-body problem can be readily simplified by 
following the approximation described in section III Bl 
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When the interactions that connects different groups 
are neglected, the evolution operator is factored out in 
groups and the spin-echo relation becomes simply 

s,w»n s ( c J> T )- ^ 

k 

where S(C g ,r) is the pseudo-spin relation, Eqn. J6J, 
for group Cg. S g (r) can be calculated numerically and 
exactly for small g (< 10). Therefore, electron- nuclear 
and intra-bath hamiltonians can be simultaneously con- 
sidered. 

In the following section, we present our algorithm 
for sorting strongly interacting nuclei in a random dis- 
tributed spin bath into well defined groupings. We 
take the electron spin-echo signal, with the initial state 
I'P) = (|0) + |l))/v2 as a figure of merit. We examine the 
convergence of our disjoint cluster approach as a function 
of the maximum group size g and consider the statistics 
of spin-echo for a variety of physical parameters such as 
Carbon-13 abundance and magnetic field magnitude and 
orientation. 



A. Grouping algorithm 

One of the criteria to aggregate groups of spins is to 
consider the strength of the intra-bath interaction. This 
parameter can be summarized in one variable C(i,j) 
which is a scalar function of the interaction C lJ between 
nuclei i and j. The aggregation algorithm used for this 
criterion is as follow. Consider an array A containing the 
criterion for all pairs ordered from high to low values in 
C and let {i,j} n be the n-th nuclear pair in array A. The 
array A is scanned completely and one of the following 
cases applies for each pair {i, j} n 

• if nuclei i & j belong to different groups: join 
both groups if TV (G(i)) + M (G (j)) < g. 

• if nucleus i belong to group G(i) and nucleus 
j does not belong to any group: add j to group 
G(z) if TV (G(z)) < g. If not, make a new group with 
j- 

• if nucleus j belong to group G(j) and nucleus 
i does not belong to any group: add i to group 
<&(j) if J\f (G(j)) < g. If not, make a new group 
with i. 

• if nuclei i & j do not belong to any group: 

make a new group with i & j, 

where Af (G) is the number of nuclei in group G and g 
is the maximum number of nuclei per group. In what 
follows, we use the criterion C (i, j) = {C l J x ) 2 + (C^) 2 + 
(C l J z ) 2 (i.e. the interaction between nuclei i and j) to 
estimate the electron spin-echo in NV centers. 
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FIG. 2: Simulation of the Pseudo spin S g =e(T) for a single 
NV center in a magnetic field of 50 Gauss oriented parallel 
the NV-axis. 
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FIG. 3: (a) Convergence: equation (|10|l as the maximum 
number of nuclei g per group is increased, (b) Indicators I g 
of the contribution of neglected pairs. When g is increased, 
the most important pair interactions are added to the pseudo 
spin relation S g . The rest is used to calculate I g . 



B. Numerical methods and example cases 

Figure [U shows S g (t) for g = 6 (Eq. ©) for 750 
random and distributed Carbon-13 in a diamond lattice 
in a magnetic field of 50 Gauss oriented along the NV- 
axis. The algorithm was implemented using MATLAB 
and the Hamiltonian for each group was diagonalized ex- 
actly followed by the calculation of the corresponding 
unitary matrices for 6000 points from to 1 ms. Each 
simulation of S g (r) takes approximately 10 min. 

The method also shows good convergence. When the 
maximum size of subgroups g is increased, more inter- 
actions among nuclei are considered and the approxima- 
tion gets better. As a figure of merit, we plot the inte- 
grated squared difference between consecutive spin-echo 
relations, S g and S g —i, 

( 5S ^ = k[ [Sg{t) - Vi(*)] a *. ( 10 ) 

Figure [3^ shows \ log (SS 2 ) up to g = 6. Each time the 
maximum size g of subgroups C g is increased, the spin 
echo relations, S g 's, get closer. 

In addition, following Ref. [TH, we introduce the follow- 
ing indicator of all interactions not included in groups 
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Cg, and therefore in S g , 



I g (r) 



{n,m}eC g 



(r) 



(11) 



The product in Eq. (|lip runs over all neglected pair inter- 
actions contained in C g and S nm is calculated according 
to Eq. ©. I g (t), which obeys < I g (r) < 1, is an 
indicator of convergence for our approach. When / s (r) 
is close to unity, good convergence is achieved. Figure 
|3b shows I g (r) for several aggregations (different g's). 
As expected, when the maximum subgroup size g is in- 
creased, the contribution from all neglected interactions 
is small. By the time the neglected interactions become 
important, the pseudo-spin S g (r) has already decayed 
(see figure [5]). 



V. RESULTS AND DISCUSSION 

The results shown in figure [4] clearly indicates that the 
electron spin echo signal cannot be modeled by just one 
time scale. This result can be understood by noting that 
few strongly interacting nuclei can coherently modulate 
the usual exponential decay. This is in good quantitative 
agreement with recent experimental results^. 

The random distribution of the spin bath and the rel- 
ative high coupling between two nearest neighbor nuclei 
(~ 2 kHz) may cause a few nuclei to contribute signifi- 
cantly to the decay of the spin-echo signal. Nuclei that 
makes small contributions to the decoherence of the elec- 
tron contribute as 1 — ar 4 exp(— ar 4 ) as it can be 
seen from Eq. ([5]). This behavior starts to deviate from 
exp(— ar 4 ) as the interaction between nuclei increases. 
Figure [^i shows a very unusual decay at which few nu- 
clei modulate coherently (Fig. rjj>, black curve) the irre- 
versible contribution from the rest of the bath (Fig. |4Jd, 
red curve). Therefore, individual NV centers can show a 
rich variety of spin-echo signals with multiple time scales. 
The coherent modulation of the spin-echo diffusion due 
to strong interacting nuclei suggests that we can think 
about a system composed of the electron and these few 
strong interacting nuclei and an environment composed 
of the rest of the spin-bath. 

Each NV center experiences a different random con- 
figuration and concentration of Carbon-13. This causes 
a large distribution of decoherence times Ti when many 
centers are probed. In order to estimate the decoherence 
time T2 we fit the envelope of S g (r) to exp (— (27-/T2) 3 ). 
When the fit is not accurate we define T2 as the longest 
time for which S g > 1/e. Figure^ shows the histogram 
of T2 for 1000 different random distribution of Carbon-13 
in the diamond lattices for an external magnetic field of 
50 Gauss parallel to the NV-axis. 

For non-zero transverse magnetic fields, the electron- 
nuclear interaction, Eq. ([7]), makes an important contri- 
bution due to the enhancement in the nuclear g-factor 
and interaction between nuclei (see Eq. (|A4|) V In this 
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FIG. 4: (a) Electron spin-echo signal highly modulated by a 
few Carbon-13 that strongly interact with the electron spin, 
(b) The strong contribution to the signal (black curve) has 
been isolated from the contribution from the rest of the spin 
bath (red curve). 
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FIG. 5: Histogram of T2 for 1000 simulations at a magnetic 
field of 50 Gauss at an angle of 6 = 0° (blue) and 9 = 6° 
(red) with respect to the NV-axis. (b) Decoherence time T2 
versus impurity concentration, Carbon-13, at 50 Gauss along 
the NV axis. 



case, nuclei will have different Larmor frequencies f2„ de- 
pending on the relative distance between them and the 
center. Each nuclei will refocus at different times, de- 
creasing the maxima of the revivals, i.e., decay of the 
signal. This is illustrated in figure EH for a magnetic field 
at an angle of = 6° from the NV-axis. 

As expected, the decoherence time decreases as the im- 
purity concentration increases. This is shown in figureEb 
where T2 goes as 1/n. To understand this it is possible 
to make an analysis using a small t expansion; while this 
is not always correct, it provides a simple explanation of 
the underlying behavior. From Eq. ([5]), the decoherence 

time scales as the geometric mean of the bath dynam- 

_ 1/2 

ics and the bath-spin interaction, i.e., T2 ~ (C^4 C ) , 
where C is the averaged nuclear-nuclear dipolar interac- 
tion and A c is some characteristic value for the electron- 
nuclear interaction. Since both interactions decay as r~ 3 
and the average distance between bodies scales with the 
concentration as n -1 / 3 , both interactions scale linearly 
in n. Therefore, the decoherence time T2 decreases ap- 
proximately as 1/n. 

As the angle between the magnetic field and the NV- 
axis, 6, is increased, the electron-nuclear dynamics dom- 
inates and the spin-echo signal shows small revivals and 
fits poorly to a single exponential decay. Thus, to de- 
scribe the coherence time at these angles, we have plot 
the average value of the signal, normalized by the average 
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f~\S B , e (t)\dt 



(12) 



Figure shows how the coherence of the signal varies 
with the strength and orientation of the magnetic field. 
This map is averaged over 6 different spin baths, since 
the random localization of Carbon-13 nuclei in the lattice 
makes the coherent time to vary from NV-center to NV- 
center as it can be seen in figure [SJs for a fixed magnetic 
field. 

When the magnetic field along the NV-axis increases, 
the contribution from the electron-nuclear interactions 
decreases (see figureJSJ:). This happens because the quan- 
tization axis for the nuclei points almost in the direc- 
tion of the external magnetic field producing a small 
oscillating field. This can be easily seen in the non- 
interacting case, Eq. (J7J), where the second term van- 
ishes if Hi "* || iln . Similarly, when electron spin-echo 
is performed using the sub manifold m s = {+1,-1}, the 
signal does not revive since each nuclei refocus the elec- 
tron at different times. This occurs because the Larmor 
frequencies in this case, fi^ 1 , are position dependent and 
differ for each nuclei. 

We also point out that the approximation introduced 
in section ["II Bl is valid as long as the impurity concentra- 
tion of Carbon-13 is not too high, so the neglected inter- 
actions that connect different groups do not play an im- 
portant role. This allows us to treat the bath as isolated 
groups. The approximation also relies on the relatively 
large interaction between the central spin (electron) and 
the bath, A n , when compared to the intra-bath interac- 
tion, C nm . The reason for this is that as the central spin 
gets disconnected from the bath (reducing A„), the de- 
cay occurs at later times r and interactions of the order 
of r _1 start to play a role. To illustrate this, consider the 
size of each subgroup scaling as {g/n) 1 ^ 3 where g is the 
size of the subgroup. Then, the interaction between near- 
est neighbor groups scales as nC nn /g where C nn is the 
nearest neighbor nuclear interaction. The time at which 
this interaction is important goes as t ~ g/nC nn = g/C. 
If we require this time to be larger than the decoherence 
time (t 3> Tq), we find that the two types of interactions 

should satisfy g [A c /Cj ^> 1 (for the present study 
this value is around 150). Therefore, when the interac- 
tion between the addressed spin and the bath is of the 
order of the intra-bath interactions, the approximation 
breaks down. This would be the case of the spin-echo 
signal for a nuclear spin proximal to the NV center- 
in which more sophisticated methods should be applied 
such as tree tensor networks^. 



VI. CONCLUSIONS 

We have presented a method to evaluate the decoher- 
ence of a single spin in the presence of an interacting 
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FIG. 6: (a) Coherent time Ti for different magnetic field 
strength and angles (measured from the NV-axis) . Each point 
is average over 6 different random distributed baths, (b) Co- 
herence time T2 versus angle of the magnetic field for 4 differ- 
ent spin baths at 50 Gauss, (c) Pseudo spin S g =6 (t) at 290 
Gauss. At high fields the collapses due to the electron-nuclear 
dynamics decreases (see text). 



randomly-distributed bath. It properly incorporates the 
strong electron-nuclear dynamics present in NV centers 
and explains how it affects the decoherence. We also 
incorporates the dynamic beyond the secular approxi- 
mation by including an enhanced nuclear g-factor that 
depend on the orientation of the external magnetic field 
relative to the NV-axis and by including an electron me- 
diated nuclei interaction. Our results show that the spin 
echo signal for NV centers can present multiple time 
scales where the exponential decay produced by many 
small nuclei contributions can be coherently modulated 
by few strongly interacting nuclei. The coherence times 
in ultra-pure diamond can be further improved by mak- 
ing isotopically pure diamond with low concentration in 
Carbon-13. This method may be used in other systems 
as long as the intra-bath interaction is smaller than the 
interaction between the central spin and the bath. These 
results have important implications, e.g., in magnetom- 
etry where long coherence times are important. For ex- 
ample, echo signals persisting for up to milliseconds can 
be used for nanoscale sensing of weak magnetic fields, as 
it was demonstrated recently^. 
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APPENDIX A: 



INTERACTIONS FOR NV 
CENTERS 



Electron spin resonance shows that the nuclear spin- 
electron spin interaction is 150 MHz 22 for the three near- 
est neighbor Carbon-13's and around 2 MHz for the Ni- 
trogen. Away from this deep defect, the interaction is 
dipolar-like, 



S ■ A n • \ n — 5.6 



^) (3(S-n)(L-n)-S-L) 



S ■ In) MHz 
(Al) 

where a nn = 1.54 A is the nearest neighbor distance for 
diamond, R n is the distance between the n-th Carbon-13 
and the defect, and n is the unit vector that connects 
the electron and the nucleus. Carbon-13's interact via 
dipolar interaction, 



this, we can write the hamiltonian in the following form 
H « AS 2 Z - le B z S z + S z A n ZJ I™ + Y, SA n- ^ 

nj n 
- 1N ]T B T • gf • I n + Y, !n • Cnm ■ lm (A3) 



where g° ff = 1 + 5g n (m s ) is the effective g-tensor 11 , 
= C nm + 5C nm (m s ) is the effective coupling be- 
tween Carbon- 1 3 's and 



T .C. .T 



2.1 



3(J„ • n)(I m ■ n) - I n ■ I m J kHz 
(A2) 



where R nm is the distance between the n-th and m-th 
Carbon-13. 

Recent experiments have been performed at fields of 
1 — 100 Gauss^i with Zeeman energies of ^MHz and 
~kHz for the electron and Carbon-13 nuclear spins, re- 
spectively. At these fields, it is convenient to make the 
so-called secular approximation thanks to the large split- 
ting A. However, transverse directions (perpendicular to 
the NV-axis) of the external magnetic field and hyperfine 
fields from Carbon-13 must be taken into account as a 
perturbation to get an accurate description of the sys- 
tem, as it has been shown by experiments^. Including 



$gn {m s 



An An An 

(2-3K|) 7e / A ™ A A Z 

yx yy yz 



A77V 



000 



SCnm (m s ) = - 



(A 7A r/ 7e ) T 

7<>gn ■ °gr, 



A(2-3|m s |) 



(A4) 



5A*(m s ) = 



(2-3K|) 
2A 



For nuclei close to the center, Sg can reach values between 
— lb 11 , SC can be several times the bare dipole-dipole 
interaction 2 ^. The term 8 A is a small contribution that 
can be neglected in most of the cases due to the large 
value of the zero-field splitting A. 
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